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Dusty disks around young stars are formed out of interstellar dust that consists 
of amorphous, submicrometre grains. Yet the grains found in comets [lj and 
meteorites [2], and traced in the spectra of young stars [3], include large crystalline 
grains that must have undergone annealing or condensation at temperatures in 
excess of 1,000 K, even though they are mixed with surrounding material that 
never experienced temperatures as high as that [4j. This prompted theories of 
large-scale mixing capable of transporting thermally altered grains from the inner, 
hot part of accretion disks to outer, colder disk regions [5], [6], [7], but all have 
assumptions that may be problematic f8ll9~1,[T0~ 1,lllU12| . Here I report that infrared 
radiation arising from the dusty disk can loft grains bigger than one micrometre 
out of the inner disk, whereupon they are pushed outwards by stellar radiation 
pressure while gliding above the disk. Grains re-enter the disk at radii where it is 
too cold to produce sufficient infrared radiation pressure support for a given grain 
size and solid density. Properties of the observed disks suggest that this process 
might be active in almost all young stellar objects and young brown dwarfs. 

^ ■ The history of thermal and compositional alternation of dust in dense dusty protoplan- 

etary disks around young pre-main-sequence (PMS) stars enables us to better understand 
conditions that initiate formation of planets. One of the long standing problems arising from 
this approach is the presence of crystalline dust in disk environments considered too cold for 
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crystallinity to occur. Thus, it has been suggested that silicates crystalize in the hot part 
of the disk close to the central star and are transported outward into colder environment. 
Currently favored theories of outward transport include: i) turbulent mixing [5], ii) ballistic 
launching of particles in a dense wind created by interaction of the accretion disk with the 
young star's magnetic field (X-wind model) [5], and iii) mixing mediated by transient spiral 
arms in marginally gravitationally unstable disks[7j. Although these theories sound promising 
and may eventually result in the definitive solution to the problem of a large scale mixing, 
they are so far hampered by theoretical assumptions needed for them to work. The turbulent 
mixing requires a source of efficient turbulent viscosity and the magnetorotational instability 
(MRI) is invoked as the most promising candidate, but large stretches of the disk are consid- 
ered not sufficiently ionized to keep MRI active [H [9j [10]. The X-wind model relies on the 
theoretical notion of magnetic field configurations in the immediate vicinity of PMS stars and 
high hopes are put on future observations to resolve this predicament [IT]. The spiral arms 
model is in the domain of discussions whether the underlying numerics, physical approxima- 
tions and assumptions on the initial conditions are realistic enough to make results plausible 

I81EIH2]. 

Unlike these theories, non-radial radiation pressure does not require additional assumptions 
on the physical conditions in the disk because it stems from the basic radiative transfer 
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properties of optically thick dusty disks. It has been already shown that individual submicron 
grains do not move far away in the disk when pushed by radiation pressure because the force 
is primarily produced by radial stellar flux [13] . On the other hand, micron or larger grains are 
large enough to also have efficient interaction with the near infrared photons (NIR; equivalent 
to dust temperatures of ~1000-2000K) from the hot inner disk. Submicron grains are very 
inefficient emitters in NIR, hence, they overheat and sublimate further away from the inner 
disk surface. This leaves the surface populated only with large grains, while small grains can 
survive within the optically thick interior [TJ] or at larger disk radii. Direct imaging with NIR 
interferometers revealed that the observed location of inner disk rim is consistent with this 
description (see [T6l [TT] and references therein). 

In optically thick protoplanetary disks dust particles < 1mm are well coupled with the 
gas and their dynamics is dominated by the gas drag [TBI [T8] . Hence, dust motion is very 
similar to the gas orbital, almost Keplerian, motion. Radiation pressure force serves as a 
slow perturbation that leads to the rearrangement of dust orbits. In order to reconstruct 
the trajectory of particles pushed by radiation we need to derive the spatial orientation of 
radiation pressure vector. For that, we need estimates of the diffuse flux as the source of 
pressure asymmetry. I solve this using the two-layer formalism, which is a well established 
method utilized in problems involving protoplanetary disk emission [19|. 

A short simplified solution is presented in figure [I] while a more rigorous derivation, which 
includes gravity, gas drag and radiation pressure, is described in Supplementary information 
§1. The result shows that the net radiation pressure force, which combines stellar and diffuse 
flux components, is directed exactly parallel to the disk surface irrespective of its curvature. 
This leads to a very interesting scenario. If the force is strong enough to move a large dust 
grain then such a large crystalline grain formed at the hot inner rim would glide over the 
disk surface toward colder disk regions until the diffuse disk flux becomes too "cold" (i.e. its 
peak wavelength is larger than the dust size), at which point the force keeping the dust afloat 
ceases. 

Further insight into the nature of non-radial radiation pressure outflow requires a more de- 
tailed description of the disk structure and an advanced radiative transfer calculation. I started 
with preliminary modeling at such an advanced level. The first results are presented in figure 
[2j The model assumes dusty disk density structure of the form pa(R,z) ~ R~ 2 exp(— z 2 /2h 2 ), 
with the scale height h = 1.67 x 10 _2 i? 1,25 , where R and z are cylindrical coordinates scaled 
with the dust sublimation radius Ri n (the disk's inner rim; see figured]). The disk contains 
0.1/zm and 2/zm olivine grains with the relative density ratio 10 4 :1 and the overall radial visual 
optical depth at z = of 10,000. I performed a full 2D radiative transfer for the case of disk 
heating from a 10,000K star and 1,500K dust sublimation. Location of the dust sublimation 
disk surface is calculated self-consistently from the mutual exchange of infrared energy be- 
tween 0.1/zm and 2/zm grains, resulting in R in = 44.7R* (R* is the stellar radius). Figure [2] 
shows the map of vertical radiation pressure on 2/im grains and examples of grain trajectories. 
Results from this detailed approach confirm the plausibility of our theoretical arguments. 

The ability of large grains to migrate along any disk curvature makes this theory indepen- 
dent of the ongoing debate on the geometrical structure of the inner disk region [15j . The 
popular view is that the inner sublimation edge is puffed up and curved [TT] . The non-radial 
pressure would affect dust dynamics under such a disk curvature in the same way as in the 
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numerical example above, except that individual grains could decouple more easily from the 
inner disk and fly toward outer disk regions due to the disk's self-shadowing [20] . 



Grains pushed by radiation create an outflow that operates at much shorter timescale than 
the local dust settling because radiation pressure is active in the region of lower gas density. 
In Supplementary information §2 I provide an estimate of the total amount of dust that flows 
outward in the disk surface layer. The outflow strength and range depend on the ratio, (3, 
of radiation force tangentially to the disk surface over local gravity force (see Supplementary 
information §1 for a detailed description): 



(3 ~ 0.4 
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where a is the dust grain radius, L* is the stellar luminosity, M* is the stellar mass and g s 
is the grain solid density. Grains with (3 > 0.5 are gravitationally decoupled from the star 
and will be pushed away from the star as long as the diffuse flux keeps them afloat within 
the optically thin surface. Grains with (3 < 0.5 feel a "reduced" gravity and their settling is 
slowed down. 



I made an attempt to estimate the spatial extent of significant vertical radiation pressure 
along the disk surface. I use a simplified, but illustrative model of the protoplanetary disk 
where the disk surface contains only single size grains. Results show (see figure [3] and Sup- 
plementary information §3) that significant dynamical effects from the non-radial radiation 
pressure are possible only for grains larger than about 1/zm. Grains a few microns in size can 
be lifted out of the disk only at small disk radii where the disk is the hottest, but already 
5/xm grains can "glide" to large radii (over 1,000 stellar radii), provided that the radiation 
pressure is strong enough to push such a grain. The upper limit on grain size pushed that 
way is dictated by equation [1] that shows how the force decreases with grain size. 

Notice that I assume a solid spherical grain, which is a simplification of a more realistic 
fluffy dust aggregate [2D [22] • Aggregates result in a much larger (3 for the same grain size 
because they have a much lower grain density than the typical 3000kg /m 3 due to inclusion of 
vacuum into the grain structure. On the other hand, crystalline grains are largely transparent 
in the spectral range of stellar radiation [lj, which would make radiation pressure ineffec- 
tive. This remains an open problem for our theory, although crystalline grains incorporated 
into dust aggregates might have a non-transparent "glue" keeping the aggregate together, 
which would increase (3 and mitigate these problems. Such "dustballs" are considered to be 
precursors of chondrules and CAIs in meteorites [23J. 

The main stellar parameter dictating the overall strength f3 of radiation pressure effect on 
a grain is the luminosity- mass ratio L*/M*. Observations and evolutionary tracks indicate 
that L^/M* > 0.5 (which gives (3 ~ 0.4 for a grain of 1/zm diameter) in almost all young stellar 
objects, including brown dwarfs. Thus, non-radial radiation pressure is at least marginally 
relevant in all these objects, especially if a realistic dust aggregate model is taken into account. 
Moreover, at earlier evolutionary stages (3 was larger because, according to stellar evolution 
models, the end of significant accretion (99% of the final mass) ends with L^/M* > 10 for 
stars M* < 1M Q [21]. 

Since crystallization is very efficient along the hot inner disk rim, radiation pressure mixing 
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of large grains would inevitably include the crystalline fraction and disperse such dust over 
the disk surface. Interestingly enough, such a correlation between large grains and crystalline 
fraction is detected in Herbig Ae stars (e.g. [25j EE]). This would be the most pronounced in 
the inner disk regions, closer to the inner rim, as it is indeed observed (e.g. [3j ETJ, EH])- With 
the help from disk turbulence, the surface of inner disk region is constantly replenished with 
new grains and the process continues as long as the radiation pressure is active. 
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Figure 1: Geometry of non-radial radiation pressure. The sketch shows a cross section of an 
optically thick protoplanetary dusty disk heated by a star. The disk has a central hole of radius Ri n 
where the dust overheats and sublimates away. According to the two-layer formalism [19], analysis 
of the disk emission in the near IR can be reduced to the disk's optically thin surface, which is 
heated directly by the stellar radiation. In this approach the disk surface is replaced with a single 
temperature layer and we assume that the stellar radiation is completely absorbed within this layer. 
The disk interior is described as the second temperature layer, but it is heated only by infrared 
radiation from the surface layer and, therefore, it is much colder and does not contribute to the disk 
emission in the near and mid IR [19] . In optically thick passive disks we can use energy conservation at 
a surface point (R, z) to set balance, F* sin a = F d , between the bolometric stellar flux F* intercepted 
by the disk at a grazing angle a and the outgoing disk radiation F d (IR emission and scattered stellar 
photons). In the approximation of geometrically thin disk surface we can assume that the entire local 
diffuse flux at the very surface is perpendicular to the surface. Grains that manage to decuple and 
move away from the surface would feel a reduced flux since the diffuse radiation streams out in all 
directions. We can decompose F d into radial, F d = — F*sin 2 a, and azimuthal, Ff_ = sin a cos a, 
components. If dust grains are big enough to have constant extinction in the wavelength range of 
F d then the radiation pressure force becomes $ cx F* + Fd- Using flux components from above gives 
the radial force 3V oc i 7 * cos 2 a and the azimuthal force $± oc F* sin a cos a. Notice that this yields 
radiation pressure force directed exactly parallel to the disk surface, $±/$r = tana, irrespective 
of the disk curvature. A more rigorous derivation is presented in Supplementary information §1, 
including dust dynamics due to gravity, gas drag and radiation pressure. 
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Figure 2: Trajectory of dust grains under the influence of stellar gravity, gas drag and 
non-radial radiation pressure. The colored background map shows the vertical z component of 
the radiation pressure vector scaled with the value that the stellar pressure would have if the dusty 
disk were not there. Two stellar luminosty-mass ratios are used: 80Lq/Mq (white dashed line) 
and 25Lq/M & (white solid line). Dust composition is olivine [29] of 2/im radius and solid density 
3000kg /m 3 . Dust grains start their travel with a vertical upward motion until the gas density drops 
enough to loosen the influence of gas drag. After that the grain is ejected to a larger disk radius, 
where trajectory details depend on the strength and direction of radiation pressure. Trajectory is 
calculated numerically with the Runge-Kutta method. Radiation pressure is calculated numerically 
from 2D radiative transfer that includes dust absorption, scattering and emission. The disk consists 
of 0.1//m and 2/xm grains that sublimate at 1500K, but the surface in this disk region is too hot 
for 0.1//m grains, which survive below the surface populated by 2/im grains. Spatial dimensions are 
scaled with the disk sublimation radius Ri n (see figure [[]) . The disk gas and dust densities decrease 
exponentially with z. Red lines show the disk surfaces defined by the radial visual optical depth 
of 0.1 (dashed red line) and 1 (solid red line). Details of this numerical result will be shown in a 
separate publication. 
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Figure 3: Estimated strength of diffuse radiation pressure along the disk surface, indi- 
cating how far grains can travel. The estimated strength, B, is denned as the ratio of diffuse 
to stellar (3 perpendicular to the disk surface (see equation 34 in Supplementary information §3 for 
details), at various distances from the star. Optical properties of the pushed grains and dust form- 
ing the disk surface are the same. The surface contains only one grain size and type. Lines show 
results for spherical grains of 0.5/um, 1/^m and 5/mi radius. Diffuse radiation pressure is important 
(i.e. B ~ 1) only for grains >l/um. Grains >5/mi experience strong diffuse pressure over a large 
disk surface because of their efficient infrared absorption at longer wavelengths, while smaller micron 
grains can float only above the inner disk with the highest temperature. The dust is enstatite from 
[30| . Other compositions lead to qualitatively similar curves. Lines start at radii defined by 1,500K 
dust sublimation temperature. 
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1 Equations for non-radial radiation pressure dynamics 



The radiation pressure force pushing a grain of radius a in direction h is 

f = —[ QTFxd\ (2) 



where c is the speed of light, Q^ 1 is the extinction coefficient and F\ is the total radiation 
flux in direction n. The stellar contribution to the radiation pressure on a grain at distance r 
from the star is 

& = ^ / QTF.xdX = ^J QTf*xd\, (3) 

where is the stellar luminosity and f*\ is normalized shape of stellar spectrum J f*\d\ = 1. 
We use 3* to scale the force 

T jQT'F.xdX jQ e x xt F*xd\- [) 

In the case of a negligible diffuse flux the pressure becomes identical to the stellar radiation 
force and f = f . 

If dust grains are big enough to have a constant extinction in the wavelength range of F x 
then the radiation pressure force becomes 5 oc F* + F d . The disk flux is perpendicular to the 
disk surface and we can decompose it into radial and azimuthal components (see figure [T]) 

F r d = — F* sin 2 a, (5) 

F'l = F* sin a cos a. (6) 
Using these components gives the radial force 

$ r oc F* — F* sin 2 a = F* cos 2 a (7) 

and the azimuthal force 

$± oc F* sin a cos a. (8) 



We see from this analysis that the radiation pressure force is directed exactly parallel to the 
disk surface $±/$ r = tana irrespective of the disk curvature. This leads to a very interesting 
scenario. If the force is strong enough to move a big dust grain then such a big crystalline 
grain formed at the hot inner rim would glide over the disk surface toward colder disk regions 
until the diffuse disk flux becomes to "cold" (i.e. its peak wavelength is larger than the dust 
size), at which point the force keeping the dust afloat ceases. 

The "strength" of the radiation pressure force is measured by the ratio of the radiation 
pressure in direction n and the local gravity force 
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which yields 



(3 = 0.2 



3000kg /m 3 

Qs 



f Qxf*\dX 



(10) 



where g s is the solid density of a grain and L & and M Q are solar luminosity and mass. Grains 
with j3 cos a > 0.5 are gravitationally decoupled from the star and will be pushed away from 
the star as long as the diffuse flux keeps them afloat within the optically thin surface. Grains 
with (3 cos a < 0.5 "feel" a reduced gravitational force and shift to a larger stable orbit. 



The pressure vector f is close to unity and, as described above, points tangentially to 
the disk surface. Stellar radiation peaks at wavelengths smaller than big grains, hence we can 
approximate Qf^ 1 ~ 2 and J Q e ^ 1 f*\d\ ~ 2. This value of Q e * 1 is correct independently of 
all other grain properties (chemical composition or shape) when the grain size is much larger 
than the wavelength (so called "extinction paradox"). In cases when the grain size is larger 
by a factor of a few, Q e ^ 1 can vary between ~ 1 and ~ 4 due to the contribution of anisotropic 
scattering. We see from this that the stellar luminosity to mass ratio and grain density dictate 
the size of a grain capable of migrating out of the hot inner disk rim. 



The equation of motion of a particle in a gaseous medium includes forces of gravity, gas 
drag in the Epstein regime and radiation pressure [131 EI] 



-G— r 



Qg_ Cs_ 

Qs a 



(11) 



where g g and v£ are the local gas density and velocity, respectively, and c s is the local speed 
of sound. We can assume that on the time scales of interest the gas has no radial or vertical 
velocities and rotates around z axis with the Keplerian angular velocity. We expand this 
equation in a cylindrical coordinate system (R, <p, z) 



R = R<p 2 - 



GM* R 

{R 2 + z 2 f/ 2 



- fiR + p 



R 



_GM 1 _ 

R 2 + z 2 ' 



GM, z 



{R 2 + z 2 fl 2 



n^R 

/iz + (3 Z 



GM* 



R 2 + z 2 



where we use = [3rR + [3 z z and 



Qg C s 

11 = -*-—. 

Qs a 



(12) 

(13) 
(14) 

(15) 



We work in the regime /i ^> R/R where particles are strongly coupled with the gas and have 
a short gas drag stopping time. Hence, in equation [T3] we can assume that dust and gas have 
the same angular motion similar to the Keplerian speed VL 2 K = GM^/R 3 pjjj. Replacing (p in 
equation [12] with Qk and using 

- 2 a m 

(16) 



GM. 



GAL R 



R? (R 2 + z 2 ) 3 / 2 



GM, 3 z?_ GM, 
R 2 2 R 2 ^^ R R 2 + z 2 



yields the solution 



R 



GM, R 
H R 2 + z 2 ' 



(17) 
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Similarly, from equation [TH we derive 



A - • (18) 



These are velocities of big dust particles in the optically thin disk surface under the influence 
of stellar and diffuse radiation pressure. 



Direction of trajectories in the R — z plane is equal to the ratio 

z (3 Z — sin 9 
R = Wr ' 

where sin 9 = zj VR 2 + z 2 . From equations [7] and [8] it follows 

(3 Z §y sin 9 + $± cos 9 sin 9 cos a + cos 9 sin a 
Pr $r cos 9 — $± sin 9 cos 9 cos a — sin 9 sin a 

In the inner disk region cos 9 cos a 3> sin # sin a, which yields 



(19) 



(20) 



z sin 

^~tan6» + tana — . (21) 

R Pr 



Notice that under a strong radiation pressure force (i.e. sm9/ ' (3r <^ 1) tan^ + tana is 
exactly the curvature of the disk surface for small angles a and 9 (figure [T|). The same is 
true for the inner disk rim where a is not small, but 9 <C 1, and the trajectory becomes 
z/R ~ tan a. 

Vertical radiation pressure f3 z reduces the influence of gravity on big grains, which results 
in expansion of optically thin disk surface. This is equivalent to a disk where gravity on a 
given grain size is reduced by a factor of 1 — /3 Z , which increases the scale height by 1/ a/1 — f3 z . 
Grains with (3 Z > 1 can decouple from the dense gaseous disk if they reach heights where gas 
drag does not dominate the dust dynamics. Such a vertical expansion works only with big 
grains and optically thin dust. // too much dust enters this zone and makes it optically thick, 
the radiation pressure decreases and the expansion subdues. 
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2 Dust mass flowing in the disk surface 



The outflow velocity is (see equations [T71 and [T8|) : 

GM* 



+ z 



-•2 







fi(z s ) R 2 + z 2 



(22) 



We assume that v is constant within the optically thin surface that starts from the height z s . 
For the surface populated by grains of average radius a the net dust mass flux is [13] 



2 / v(R,z)g d (R,z) 2nRdz, 



(23) 



where Qd(R,z) is the dust number density at (R,z). The factor of two comes from the disk 
having two sides. Using // from equation [15] and gas to dust ratio £ = g g / g^, the mass flux 
becomes 

oo 

GM*Q s a(3 „ f dz 



47T- 



-R 



R 2 + z 2 



(24) 



The assumption of a constat sound speed c s ~ 2000m/s is correct within about ±700m/s for 
temperatures considered here. Solving the integral and neglecting z s /R <C n/2 yields 



0.021 Ma 



£ year 

We replace (3 with equation [10] and get 

g 



M, 



3000A; g/m 3 



jj,m 



(3. 



x 10" 3 M ff 



i year 

where we assume J Q^f^xdX ~ 2 and | f | ~ 1. 



(25) 



(26) 



Estimated amount of dust from equation [26] that flows outward in the disk surface layer 
depends critically on the gas to dust ratio £, which is very uncertain in the inner disk region 
because of the interplay between dust sublimation, growth and settling. For the standard 
£ = 100 the outflow transfers one Earth mass within ~13,000 years for = 1L Q and ~13 
years for L* = 100L Q . If the dust outflow becomes larger than the disk accretion inflow then 
the dusty disk starts to erode from its inner rim outward. According to our estimation, such 
an erosion happens when the total disk accretion is < 2.4 x lO _8 [L^/L ]M Q /yr. This value 
is about the same as the observed accretion rate averages in T Tauri [32] and Herbig Ae [33] 
stars, while in more luminous Herbig Be stars the limit becomes very high. Interestingly, 
observations indicate structural differences in the inner disk geometry between low and high 
luminosity young stellar objects [TB] [IB"]. 

The timescale for an outflow of dust from cylindrical radius R to R' > R can be derived 
from equation [22] The radial component of the velocity is 



GM.P 



R 



GM,(3 J_ 

H(z s ) (R 2 + z 2 )V 2 ~ y,(z s ) R 2 ' 



(27) 
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which gives the time for a grain transport from RtoR'>R 



R' 



trad(R, R) 



dR/v R 



3GM*(3 



R 



R*_ 
~R 



fJ>(z s 



3pn 2 K (R) 



R* 



- 1 



(21 



where fj,(z) is defined in equation [151 z s is the height where optically thin surface starts and 
VLk{R) is the Keplerian speed at R. In comparison, the local dust settling time is t se t{R) = 
li/VL 2 K (R) [IB], where \i goes over all z, yielding t set ~ 10 5 yr for \im sized particles in typical 
disks. Radiation pressure is active in the region z ~ z s where /i is much smaller than in 
the disk interior of small z. Therefore, it operates at much smaller timescales than t set when 
R' ~ R, but becomes comparable t rad (R, R') ~ t set when R'/R > 10. 
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3 Extent of the diffuse radiation pressure 



We compare values of j3 perpendicular to the disk surface originating from diffuse {(3ir) and 
stellar radiation (/3* sin a) 

(3* sin a J Q' x F*d\ sin a ' v ; 

where Q' x is the extinction coefficient of a dust grain pushed by the radiation pressure. For a 
given distance r from the star, the stellar component can be rewritten as 

Q' x F* x d\=(^j a SB T*\Q') T », (30) 

where <jsb is the Stefan-Boltzmann constant, R* and T* are the stellar radius and temperature, 
respectively, and {Q')t* is averaged Q' x over the stellar spectrum. 

The diffuse radiation is difficult to specify because it depends on structural properties of a 
particular protoplanetary disk, hence it suffers from various modeling uncertainties. But we 
know that the most inner dusty disk region contains only big grains in its hot surface layer 
because small grains overheat at these distances and sublimate under direct stellar radiation. 
Hence, here we consider a disk surface populated with big grains. The radial thickness of disk 
surface is defined by complete absorption of stellar radiation, where we use the surface radial 
optical depth at A=0.55/zm as ry ~ 1. The diffuse radiation is approximated according to the 
approach described in [19] 



/ Q' x F x d\ = a SB T A J Q x b x (T) e IR sma d\, (31) 

where T is the surface dust temperature, b x (T) is normalized Planck function and e/^sino; 
is the surface thickness in infrared. Unlike small grains where em = Q^/Qy 8 < 1 always 
holds, big grains have 

Qf s /Qv s when Qf s < Of 



IR \ 1 when Qf s > Q^ s 1 ; 

where Qf"* is the absorption coefficient of surface dust. The surface temperature is dominated 
by stellar heating, hence we can estimate the temperature from the equilibrium between stellar 
absorption and optically thin infrared emissioiJl] 

T\ A (R*\ 2 (Q abs ) T * 



T* J \r J (Q abs ) T ' (33) 
where (Q abs )T is the Planck average at temperature T. Combining 1301 [3X1 and [331 with [291 gives 

B = Pir = {Q abs )T* (Q'e IR ) T , . 

(3* sin a (Q abs ) T {Q'}t* ' 1 ' 



1 Equation [33] holds for the optically thin surface of an optically thick disk, while optically thin disks have 
a factor of 1/4 in the right hand side of the equation. For the source of this difference between optically thin 
and thick disks see [21 [51] ■ 
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